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The  Meteorology  Division,  Air  Force  Geophysics  Laboratory,  operates 
an  MC130E  cloud  physics  research  airplane  which  is  outfitted  with  a variety 
of  instruments  to  measure  hydrometeor  concentrations  aloft.  An  imoortant 
application  of  this  research  airplane  is  to  determine  concentrations  of 
hydrometeors  which  could  cause  missile  erosion. 

The  instrument  of  concern  here  is  called  the  EWER  (Evaporation  of  Water 
that  Erodes  on  Reentry)  which  measures  hydrometeor  water  content.  It  is 
provided  with  a relatively  large  intake  aperature  (10  cm2),  so  to  adequately 
sample  low  concentrations  of  hydrometeors,  which  is  analogous  to  the  wide 
mouth  of  a ewer  jar. 

Calculated  collection  and  measurement  efficiencies  of  the  EWER  are  re- 
ported here  for  water  drops  and  five  types  of  ice  hydrometeors  over  wide 
ranges  of  particle  sizes  and  flight  conditions. 
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EWER  DESCRIPTION 


In  external  appearance  the  EWER  consists  of  a pair  of  tapered  tubes 
(probes)  protruding  from  the  forward  end  of  the  left-hand  wing  pod  on  the 
MC130E  airplane  (Fig.  1).  Details  of  its  development,  construction  and 
calibration  are  given  by  Durran,  et  al.^  Figure  2 shows  critical 
dimensions. 

A schematic  which  illustrates  operation  of  the  instrument  is  "shown 
in  Fig.  3.  In  the  reference  probe  hydrometeors  are  separated  from  the 
flow  such  that  only  moist  air  passes  through  to  the  detector.  Air  and 
hydrometeors  are  allowed  to  pass  through  the  main  probe,  but  they  are 
heated  such  as  to  vaporize  all  water.  Water  content  (mass  of  water  per 
unit  volume  of  air)  is  derived  from  light  attenuation  by  the  water  vapor 
which  passes  through  the  detector.  The  detector  consists  of  four  Lyman-o 
humidometers  mounted  on  a rotating  head  such  that  each  humidometer  in  turn 
samples  both  the  total  and  reference  water  vapor.  The  reference  water 
vapor  content  is  subtracted  from  the  total  water  vapor  content  to  yield 
the  hydrometeor  water  content. 

Air  sampling  by  the  EWER  is  not  isokinetic.  Air  sampling  efficiency, 
n (airflux  into  the  EWER  divided  by  free  stream  air  flux),  is  shown  in 
Fig.  4 as  a function  of  flight  (i.e. , free  stream)  airspeed.  For  our 
calculations  we  have  taken  n to  be  0.45. 


1.  D.  A.  Durran,  D.  H.  Ross  and  W.  J.  Swartwood,  "Development  of  Cloud 
Water  Content  Meter,  EWER,"  The  Aerospace  Corp. , Report  SAMSO-TR-78- 
113  (March  1978).  AD-A061  255 
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Figure  1.  Wing  pod  and  EWER  mounted  under  the  left  wing  on  the  MC130E. 


Figure  3.  EWER  Schematic.  (From  Durran,  et  al.'  ') 


COLLECTION  AND  MEASUREMENT  EFFICIENCY  CALCULATIONS 


GENERAL  DISCUSSION 


The  EWER  probes  are  mounted  far  eivgh  from  the  airplane  wing, 
propellers  and  fuselage  that  their  orifices  are  essentially  in  the  free 
stream.  However,  as  noted  above  air  sampling  is  at  about  45%  of  the  free 
stream  Hux  under  normal  flight  conditions.  Since  partial  stagnation  occurs 
at  the  orifices,  the  question  arises  as  to  how  this  affects  hydrometeor 
sampling  and  measurement. 

To  calculate  collection  efficiency,  we  develop  a digital  description 
of  the  EWER  probes  and  wing  pod  nose  and  use  this  to  calculate  air  flow 
around  and  into  the  EWER  orifices  for  specified  flight  conditions.  Then  we 
compute  trajectories  of  hydrometeors  through  the  three-dimensional  flow 
field  from  the  undisturbed  cloud  into  the  orifices  such  as  to  define  a 
quantity  called  concentration  factor. 

Concentration  factor,  Cp,  is  defined  as  the  ratio  of  particle  flux  at 
the  sampling  or  target  point,  Ft,  to  the  particle  flux  in  the  free  stream,  F, 


C 


F = -r  . 


(1) 


Concentration  ratio,  CM  t,  is  the  ratio  of  the  particle  concentration  at  the 
target  point  to  free  stream  concentration. 


CM.t  * CF  Wp.t  > <2> 


where  V . is  the  particle  speed  into  the  orifice  and  V is  free  stream 

P » L 

airspeed. 

2.  H.  G.  Norment  and  R.  G.  Zalosh,  "Effects  of  Airplane  Flowfields  on  Hy- 
drometeor Concentration  Measurements,"  AFCRL-TR-74-0602  (6  Dec.  1974) 
AD-A006  690 

3.  H.  G.  Norment,  "Effects  of  Airplane  Flowfields  on  Cloud  Water  Content 
Measurements,"  AFCRL-TR-75-0231  (30  April  1975).  AD-A014  807 

4.  H.  G.  Norment,  "Additional  Studies  of  the  Effects  of  Airplane  Flowfields 
on  Hydrometeor  Concentration  Measurements,"  AFGL-TR-76-0187  (13  August 
1976).  AD-A032  311 
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Concentration  factor  is  determined  by  computing  particle  trajectories 
from  the  free  stream  (initial  plane  in  Fig.  5)  to  a small  area  in  the  target 
plane  (Fig.  5)  that  surrounds  the  target  point  such  as  to  define  a particle 
flux  tube.  Since  the  particle  mass  transfer  rate  through  the  tube  is  con- 
stant at  all  cross  sections,  it  is  easily  shown  that 


where  A and  At  are  the  flux  tube  cross  sectional  aceaS'in'fche  freC-stream  I 
and  at  the  target  point.  In  the  limit  as  A and  At  approach  zero,  eq.  (3) 
becomes  exact. 

For  the  situation  considered  here  the  target  plane  is  constrained  to  lie 

parallel  to  the  EWER  orifice  plane  at  a distance  1.5  mm  upstream  from  it. 

Target  points  are  defined  near  the  orifice  periphery,  at  95%  of  the  orifice 

radius,  such  that  the  flux  tube  incident  on  the  target  plane  envelopes 

approximately  90%  of  the  flux  into  the  orifice. 

The  water  content  implied  by  eq.  (2)  applies  only  to  the  EWER  orifices; 

it  is  not  the  water  content  seen  by  the  Lyman-a  detector  which  is  located 

inside  the  wing  pod.  This  is  because  air  entering  a EWER  orifice  has  speed 

nV  (n  = 0.45  in  this  case)  while  hydrometeors  enter  at  speed  V . > nV. 

p,t 

Once  inside  the  EWER  tube,  however,  the  hydrometeors  are  vaporized  and 
henceforth  assume  the  same  speed  as  the  air.  If  we  ignore  compression 
effects  and  recognize  that  pressure  and  temperature  differences  between 
the  detector  chamber  and  ambient  are  corrected  for  elsewhere,  a simple 
mass  balance  shows  that  the  concentration  ratio  in  the  detector  chamber, 

CM,E  1S 


= Cp/n  . 


(4) 


Thus,  while  CM  t (eq.  (2))  indicates  hydrometeor  collection  efficiency 
by  the  EWER,  the  hydrometeor  measurement  efficiency  is  given  by  eq.  (4). 
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Figure  5.  Perspective  view  of  a hydrometeor  flux  tube.  The  initial 
and  target  planes  are  perpendicular  to  the  central  tra- 
jectory (dashed  line).  Since  the  peripheral  trajectories 
are  constrained  to  pass  through  a circle  about  the  central 
trajectory  in  the  target  plane,  in  general  they  pass  through 
a non-circular,  closed  curve  in  the  initial  plane. 
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PARTICLE  TRAJECTORY  CALCULATION 

We  assume  that  the  bulk  air  flow  is  not  perturbed  by  the  particles. 
Moreover,  since  particle  density  is  large  compared  to  that  of  air,  we  can 
neglect  buoyancy  and  inertial  reaction  of  the  fluid  to  obtain  the  three- 
dimensional,  normalized  equation 


^ = f if  <VVrV-*l 

dx  fn  )vs  a p BN,sRN  ' 


(5) 


Non- dimensioned  quantities  are: 


V Va 


fn  = V2/ (Lg) 


particle  and  air  velocities 

still-air  terminal  settling  speed 

of  the  particle 

unit  vector  in  the  z (upward) 

direction 

time 

Froude  Number 


Rn  = ^ v 4 I v 

N n a p 


Reynolds  Number 


Dimensioned  quantities  are: 


$ 


! 


6 

P 

n 

g 

v 

L 


particle  dimension 
air  density 
air  viscosity 

gravity  acceleration  constant 
freestream  airspeed 
a characteristic  dimension  of  the 
fuselage 


Here  length  is  normalized  by  L,  velocity  by  V and  time  by  L/V.  RN 

N*  S 

and  Bn  s are  for  still-air,  terminal  settling  of  the  particles. 

Starting  at  the  initial  plane,  eq.  (5)  is  integrated  with  re- 
spect to  time  in  three  dimensional  space  via  the  code  DVDQ  of  Krogh^ 
until  the  target  is  reached.  The  method  used  to  compute  va  at  each  time 

a 

step  is  described  next,  and  then  drag  coefficients  are  discussed. 


THREE-DIMENSIONAL  FLOW  CALCULATION 


In  performing  concentration  factor  calculations  for  sam- 
pling sites  on  particular  airplanes,  it  is  important  to  use  three- 
dimensional  airflow.  This  is  the  only  way  to  properly  account  for: 
airplane  geometry  and  angle-of-attack,  airspeed,  altitude,  and  par- 
ticle settling. 

Cloud  physics  airplanes  are  subsonic,  sampling  runs  being 
made  typically  between  100-150  kts  indicated  airspeed  (IAS).  Par- 
ticle measurement  points  are  beyond  the  skin- friction  boundary  layer, 
and  should  be  placed  to  avoid  regions  of  separated  flow.  Therefore, 


5.  F.T.  Krogh,  "Variable  Order  Integrators  for  Numerical  Solutions  of 
Ordinary  Differential  Equations,"  Jet  Propulsion  Lab  Technology 
Utilization  Document  No.  CP-2308  (November  1970). 
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potential  (i.e.,  frictionless,  incompressible,  laminar)  flow  calcula- 
tions are  quite  adequate.  We  use  a code  developed  by  Hess  and  Smith^’^ 
for  calculating  potential  flow  about  arbitrary  three-dimensional  bodies. 
The  Hess-Smith  code  requires  input  of  a digital  description  of  the  body 
surface.  This  consists  of  the  coordinates  of  the  corner  points  of  a 
large  number  of  contiguous,  plane,  quadrilaterals.  Figure  6 shows  a 
computer  plot  of  the  digital  description  of  the  EWER  probes  and  wing 
pod  nose  used  for  this  study. 

Accuracy  of  the  flow  calculations  has  been  checked  with 
excellent  results  by  Hess  and  Smith^  for  many  bodies  for  which  ana- 
lytical solutions  are  available.  Norment  and  Zalosh^  have  compared 
computed  trajectories  around  ellipsoids  in  analytical  flow  fields  with 
similar  trajectories  in  Hess-Smith  flow  fields,  and  they  have  compared 
current  trajectory  results  with  prior  work;  agreement  is  excellent. 


PARTICLE  DRAG  COEFFICIENT 


Davies^shows  that  still-air  terminal  settling  of  spheres 
can  be  generalized  in  terms  of  the  dimensionless  numbers  RN  and 
. Over  the  range  from  the  smallest  spheres,  which  settle  under 
viscous  flow  conditions  and  obey  Stokes  law,  to  spheres  much  larger 
than  of  interest  here,  and  for  any  Newtonian  fluid,  a reproducible 
single-valued  relationship  between  R^  and  exists.  Furthermore, 
Bn  is  independent  of  settling  speed,  being  a function  of  fluid  and 
sphere  properties  only;  thus  for  given  sphere  and  fluid,  Rw  and 
hence  V$  can  be  calculated.  Polynomials  by  which  R^  can  be  com- 
puted as  a function  of  BN  $ were  derived  by  Davies  from  a composite 
of  many  sets  of  experimental  data. 


6.  J.L.  Hess  and  A.M.O.  Smith,  "Calculation  of  Non-Lifting  Potential 
Flow  About  Arbitrary  Three-Dimensional  Bodies,"  McDonnel  Douglas 
Report  E.S.  40622  (15  March  1962).  AD-282  255. 

7.  J.L.  Hess  and  A.M.O.  Smith,  "Calculation  of  Potential  Flow  About 
Arbitrary  Bodies,"  Progress  in  Aeronautical  Sciences,  Vol.  8,  edited 
by  D.  Kuchemann  (Pergammon  Press,  New,. York,  1967). 

8.  C.N.  Davies,  "Definitive  Equations  for  the  Fluid  Resistance  of 
Spheres,"  Proc.  Phys.  Soc.  (London)  57,  259-270  (1945). 
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Since  the  work  of  Davies,  it  has  been  found  repeatedly  that 
this  treatment  is  applicable  to  particles  of  other  shapes,  providing 
settling  is  steady  and  particle  orientation  is  stable. 

For  the  trajectory  calculations  required  here,  the  problem 
must  be  turned  around.  In  addition  to  gravity  settling,  there  is  a 
particle  velocity  component  (relative  to  air)  caused  by  the  distur- 
bance of  the  passing  airplane.  At  any  time  step  in  the  numerical 
integration  of  eq.  (5),  va-vp  (and  hence  R^)  is  known,  and  must  be 
determined.  For  viscous  motion  (i.e. , Stokes  flow,  where  < 1) 

= 24  Rn  and  eq.  (5)  can  be  integrated  without  question.  However, 
for  larger  R^  the  steady-state  drag  data  determined  experimentally  for 

terminal  settling  must  be  used  to  compute  accelerative  particle  motion. 

(Ql 

Experimental  measurements  by  Keimv  ' and  a theoretical  analy- 
sis by  Crowe,  et  al . (^indicate  that  if  the  acceleration  modulus. 


is  smaller  than  about  10-2,  steady-state  drag  coefficients  can  be  used 


without  significant  error  to  compute  accelerative  motion.  has 

never  been  found  to  exceed  10-2  in  our  trajectory  calculations. 

For  small  water  drops,  which  are  spherical,  the  polynomial 

equations  of  Davies^  are  used  to  compute  V.  while  for  larger,  dis- 

s (ill 

torted  drops,  the  equations  of  Foote  and  du  Toitv  ' are  used.  To  com- 


pute Bw  from  Rw,  inverse  polynomials  have  been  developed,  using  the 

M N (121 
data  set  given  by  Davies  for  small  drops  and  Gunn  and  Kinser'  ' for 


large  drops. 


9.  S.R.  Keim,  "Fluid  Resistance  to  Cylinders  in  Accelerated  Motion," 

J.  Hydraulics  Div. , Proc.  Amer.  Soc.  Civil  Eng.,  6,  paper  1113  (1956). 

10.  C.T.  Crowe,  J.A.Nicholls  and  R.B.  Morrison,  "Drag  Coefficients  of 
Inert  and  Burning  Particles  Accelerating  in  Gas  Streams,"  Ninth  Symp. 
(Int'l.)  on  Combustion,  Academic  Press,  pp.  395-405  (1963). 

11.  G.B.  Foote  and  P.S.  du  Toit,  "Terminal  Velocity  of  Raindrops  Aloft," 
J.  Appl . Meteor.  8,  249-253  (1969). 

12.  R.  Gunn  and  G.D.  Kinser,  "The  Terminal  Velocity  of  Fall  for  Water 
Droplets  in  Stagnant  Air,"  J.  Meteor.  6,  243-248  (1949). 
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Concentration  factors  also  are  calculated  for  five  different 
ice  crystal  forms  and  for  broad  ranges  of  ice  crystal  sizes.  Tables 
of  ice  crystal  properties  used  in  this  study  are  presented  below  in 
Tables  2 through  5. 
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EWER  DIGITAL  DESCRIPTION  AND  FLUX  TUBE  GEOMETRY 

The  digital  description  shown  in  Fig.  6 was  developed  from  half-scale 
engineering  drawings.  Coordinate  measurements  are  accurate  to  about  one 
millimeter.  As  shown  in  Fig.  7 the  EWER  orifice  plane  is  defined  by  twenty- 
one  quadrilaterals.  Each  of  these  quadrilaterals  is  required  to  leak  inward 
45%  of  the  free  stream  flow.  Flow  into  and  around  the  orifice  is  shown  in 
Fig.  8. 

The  target  plane  used  to  calculate  concentration  factors  is  parallel  to 
the  orifice  plane  and  offset  from  it  in  the  downstream  direction  by  1.5  mm. 
The  purpose  of  the  offset  is  to  avoid  direct  contact  with  quadrilateral  edge, 
corner  and  gap  discontinuities. 

Hydrometeor  flux  tubes  are  defined  by  seven  trajectories  for  all  cases: 

one  trajectory  to  the  center  of  the  orifice  plus  six  equally  spaced  around 

its  periphery  at  distances  from  the  center  of  95%  of  the  orifice  radius. 

Trajectories  to  these  specified  points  were  calculated  by  the  iterative 

(2) 

method  described  by  Norment  and  Zalosh.  ' Figures  9 and  10  illustrate 
trajectory  flux  tubes  for  large  and  small  hydrometeors  Respectively. 

FLIGHT  CONDITIONS 

Indicated  airspeed  of  150  knots  with  4°  angle-of-attack,*  nose  up, 
were  assumed  for  all  cases,  and  calculations  were  done  for  three  altitudes: 

2 km,  5 km,  and  7.5  km  relative  to  sea  level.  Atmospheric  conditions  at 
these  altitudes  were  taken  from  the  U.S.  Standard  Atmosphere,  1962^13^  and 
are  presented  in  Table  1 along  with  true  airspeeds. 

HYDROMETEOR  PROPERTIES 

Selection  and  analysis  of  hydrometeor  drag  and  dimension  data  are 
discussed  in  detail  in  refs.  2 and  3.  Ice  hydrometeor  properties  are 
presented  here  in  Tables  2 through  5.  In  the  trajectory  calculations  the 

-jf — — 

The  4 angle-of-attack  is  somewhat  arbitrary  in  that  the  actual  angle 
varies  slightly  with  altitude,  load  weight  and  other  flight  conditions, 
but  small  changes  in  the  angle  would  be  expected  tonhave  minor  effects 
on  the  results. 
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PROPERTIES  OF  COLUMNAR  ICE  CRYSTALS 
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TABLE  3 


PROPERTIES  OF  HEXAGONAL  PLATE  (Pla)  ICE  CRYSTALS 


Diameter  of 


Diameter 

(urn) 

Thickness 

_ipm) 

Mass 

(pg) 

Mater  Drop  of 
Equal  Mass 
(pm) 

Density 

(kg/m3) 

10 

5.68 

-v  2.582  x 10“4 

v 7.9 

* 700. 

30 

9.30 

* 3.591  x IO-3 

•v  19 

•v  660.6 

50 

11.70 

-v.  9.203  x 10-3 

% 26 

* 484.0 

100 

15.97 

* 3.351  x 10"2 

40 

* 323.0 

300 

26.16 

0.4779 

97 

312.5 

500 

32.90 

2.145 

160 

401.5 

1000 

44.91 

20.22 

338 

693.1 

- % § • 

U 


i. 

TABLE  4 

PRONWrtM  OF  PLANE  DENDRITIC  (PJe)  ICEJWSTAtS 

..  '*') 


Dendrite 

Diameter 

(pm)_ 

TMdmess 

(ym) 

Mass 

(yR) 

Diameter  Of 
Water  Drop 

Of  EfHttalMt 
Mass  (urn) 

Equivalent 

7 « * “ 

Density 

(kg/rn?) 

iUt£ 

TioTT 

Diameter 

(ym) 

500 

29.16 

'0/3366 

...  .112  - • 

M 

&1;d 

600 

31.24 

T.OMD 

128 

357 

389.4 

800 

34.82 

2.1850 

161 

402 

490.3 

1000 

37.87 

3.6484 

191 

450 

574.1 

1500 

44.13 

8.2886 

251 

556 

720.8 

1800 

47.27 

11.494 

280 

607 

785.4 

2500 

53.50 

* 21.129 

~ 343 

~ 712 

924.1 

TABLE  5 


PROPERTIES  OF  AGGREGATES  OF  UNRIMED  RADIATING  ASSEMBLAGES  OF 
PLATES,  SIDE  PLANES,  BULLETS  AND  COLUMNS 
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menslon 

Mass 

Density 

(kq/m3) 

Diameter  of 
Water  Drop  of 
Equivalent  Mass 

(ym) 

(*g) 

3.756 

(ym) 

300 

265.7 

193 

500 

9.914 

151.5 

266 

800 

24.21 

90.3 

359 

1000 

37.0 

70.7 

413 

2000 

138.1 

33.0 

641 

3000 

298.3 

21.1 

829 

. . » Jl 


hydrometeors  are  allowed  to  settle  at  their  terminal  gravity  speeds. 

We  assume  that  columns  and  plates  always  orient  such  that  the  drag 
vector  Is  normal  to  the  long  axes  of  the  columns  and  normal  to  the  faces 
of  the  plates.  These  are  considered  to  be  good  assumptions  for  this 
study. 
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RESULTS 


Results  are  presented  in  Tatles  6-10  and  CM  E vs.  hydrometeor  mass 
is  plotted  in  Fig.  11.  Comparison  of  Figs.  9 and  10  against  Fig.  8 illus- 
trates that  the  smaller  particles  tend  to  follow  the  air  flow  whereas  the 
larger  particles  substantially  ignore  the  air  flow.  Thus,  while  the  larger 
particles  enter  the  EWER  orifices  at  nearly  free  stream  speeds  and  at 
essentially  free  stream  concentrations  (see  the  V„  */V  and  Cu  * columns  in 

P rl  t C 

Tables  6 - 10),  after  vaporization  their  concentrations  approach  the 
limiting  value  of  free  stream  concentration  divided  by  n ( i . e. » the' heavy 
particle  asymptote  in  Fig.  11).  Since  the  smaller  particles  tend  to  follow 
the  air  flow,  their  concentrations  at  the  orifices  also  deviate  only 
slightly  from  free  stream  values.  However,  since  they  enter  the  orifices 
at  substantially  less  than  free  stream  speeds,  their  concentrations  after 
vaporization  also  are  less  as  shown  by  their  CM  E values,  which  in  the  limit 

of  infinitesimal  uarticle  mass  should  approach  unity. 

Durran,  et  al.^  calculated  collection  efficiencies  for  several  sizes 
of  water  drops  for  several  air  flux  efficiencies  for  atmospheric  conditions 
at  6.1  km  altitude.  Their  case  that  matches  one  of  ours  closest  is  for 
20  pm  diameter  drops  with  n = 0.6.  They  present  graphs  of  particle  flux 
ancf  concentration  across  the  orifice;  values  given” here  are  approximate- 
averages : 


Durran  et  al.^  ^ 1.06  0.85 

This  study  1.04  0.84 

Considering  that  our  calculation  was  done  for  5 km  altitude  and  for  n * 0.45, 
the  agreement  is  probably  good.  
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TABLE  6 


fl 


WATER  DROP  CONCENTRATION  RATIOS  IN  THE  EWER 


2 km 


Results  are  given  for  three  altitudes  In  the  sequence:  5 km 

7.5  km 


Water  Drop 

Diameter 

Mass 

W*  CM.t 

r 

(Mm) 

(wg) 

LM,E 

. '.**a~* 

0.626  1.070 

1.489 

10 

5.236x10-** 

0.655  1.068 

1.554 

0.682  1.064 

1.613 

- i * » 

V » 

0.777  1.048 

1.808 

20 

4. 189xl0"3 

0.804  1.041 

1.861 

0,827  1.036 

1.904 

0.847  1.035 

1.947 

30 

.01414 

0.869  1.029 

1.987 

0.886  1.024 

2.018 

0.909  1.023 

2.066 

50 

.06545 

0.924  1.018 

2.089 

0.936  1.013 

2.107 

0.945  1.012 

2.126 

80 

.2681 

0.954  1.009 

2.140 

0.961  1.006 

2.150 

0.957  1.009 

2.144 

100 

.5236 

0.964  1.006 

2.155 

0.969  1.004 

2.162 

0.983  1.000 

2.184 

300 

14.14 

0.985  0.999 

2.187 

0.987  0.998 

2.189 

0.987  0.998 

2.191 

500 

65.45 

0.988  0.998 

2.192 

0.989  0.998 

2.193 

0.990  0.998 

2.195 

800 

268.1 

0.990  0.998 

2.196 

0.991  0.997 

2.196 

0.991  0.997 

2.196 

1000 

523.6 

0.991  0.997 

2.197 

0.992  0.997 

2.197 

L. 


3 


Vp  t/V  here  and  In  Tables  7 - 10  Is  calculated  at  a point  1.5  mn  upstream  of 

the  orifice  center.  The  airspeed  ratio  at  this  point  Is  about  0.49  compared 
with  0.45  In  the  orifice  plane.  _ 
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TABLE  7 


ICE  COLUW  CONCENTRATION  RATIOS  IN  THE  EWER 

2 km 

Results  are  given  for  three  altitudes  In  the  sequence:  5 km 

7.5  km 


Column 

Solid  Columns 

Hollow  Columns 

Length* 

V t/v 

p *t 

CM,t 

CM,E 

V JM 

P»t 

CM,t 

CMfE 

0.644 

1.072 

1.535 

0.576 

1.056 

1.353 

20 

0.670 

1.071 

1.595 

0.594 

1.062 

1.402 

0.694 

1.067 

1.647 

0.611 

1.066 

1.447 

0.723 

1.077 

1.729 

0.642 

1.087 

1.550 

30 

0.753 

1.068 

1.788 

0.670 

1.083 

1.612 

0.780 

1.060 

1.836 

0.695 

1.078 

1.666 

0.839 

1.050 

1.957 

0.768 

1.073 

1.830 

50 

0.864 

1.041 

1.998 

0.799 

1.061 

1.884 

0.883 

1.034 

2.030 

0.824 

1.052 

1.928 

0.908 

1.030 

2.079 

0.858 

1.049 

2.001 

100 

0.925 

1.023 

2.103 

0.882 

1.040 

2.037 

0.938 

1.018 

2.120 

0.900 

1.032 

2.065 

0.946 

1.016 

2.137 

0.914 

1.030 

2.092 

300 

0.957 

1.011 

2.150 

0.930 

1.023 

2.114 

0.964 

1.008 

2.160 

0.942 

1.018 

2.130 

0.957 

1.012 

2.152 

0.930 

1.024 

2.116 

500 

0.965 

1.008 

2.162 

0.943 

1.018 

2.134 

0.971 

1.005 

2.169 

0.953 

1.013 

2.146 

0.968 

1.007 

2.166 

0.948 

1.016 

2.142 

1000 

0.974 

1.004 

2.173 

0.958 

1.011 

2.154 

0.978 

1.002 

2.178 

0.965 

1.008 

2.163 

0.978 

1.002 

2.179 

0.966 

1.008 

2.165 

3000 

0.982 

1.001 

2.183 

0.972 

1.005 

2.172 

0.984 

1.000 

2.186 

0.977 

1.003 

2.177 

♦ 

See  Table  2 for  ice  column  properties. 
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TABLE  8 


HEXAGONAL  PLATE  (Pla)  CONCENTRATION  RATIOS  IN  THE  EWER 


2 km 


Results  are  given 

for  three  altitudes  In  the  sequence: 

5 km 
7.5  knr 

Plate 

Diameter* 

(um) 

Vt/V 

CM.t 

CM,E 

0.643 

1.062 

1.518 

10 

0.676 

1.060 

1.594 

0.708 

1.056 

1.661 

0.718 

1.075 

1.715 

30 

0.750 

1.067 

1.779 

0.779 

1.059 

1.833 

0.733 

1.078 

1.756 

50 

0.764 

1.069 

1.817 

0.792 

1.061 

1.867 

0.753 

1.083 

1.812 

100 

0.784 

1.073 

1.869 

0.811 

1.063 

1.916 

0.842 

1.062 

1.986 

300 

0.869 

1.060 

2.027 

0.890 

1.040 

2.058 

0.890 

1.043 

2.064 

500 

0.911 

1.033 

2.092. 

0.926 

1.026 

2.112 

0.945 

1.019 

2.140 

1000 

0.956 

1.013 

2.153 

0.964 

1.009 

2.162 

See  Table  3 for  properties  of  hexagonal  plates. 
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TABLE  9 


PUNE  DENORITE  (Pit)  CONCENTRATION  RATIOS  IN  THE  EWER 

2 km 

Results  are  given  for  three  altitudes  In  the  sequence:  5 km 

7.5  km 


Dendrite* 

Diameter 

(pm) 

V* 

CM,t 

CM,E 

0.861 

1.053 

2.016 

500 

0.886 

1.042 

22052 

0.904 

1.034 

1.079 

0.895 

1.040 

2.070 

800 

0.915 

1.031 

2.096 

0.930 

1.024 

2.116 

0.910 

1.035 

2.092 

1000 

0.927 

1.026 

2.114 

0.940 

1.020 

2.131 

0.940 

1.020 

2.133 

1800 

0.952 

1.015 

2.148 

0.961 

1.011 

2.158 

0.952 

1.015 

2.148 

2500 

0.962 

1.010 

2.160 

0.968 

1.007 

2.167 

*See  Table  4 for  plane  dendrite  properties. 
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TABLE  10 

CRYSTAL  AGGREGATE  CONCENTRATION  RATIOS  IN  THE  EWER 


Aggregate 

Dimension* 

J*L- 


for  three  altitudes 

In  the  sequence: 

5 km 
7.5  km 

vp.t/v 

CM,t 

CM,E 

0.970 

1.003 

2.162 

0.974 

1.001 

2.169 

0.978 

1.000 

2.174 

0.975 

1.001 

2.170 

0.978 

1.000 

2.175 

0.981 

0.999 

2.179 

0.978 

1.000 

2.175 

0.981 

0.999 

2.179 

0.983 

0.999 

2.182 

0.980 

1.000 

2.177 

0.982 

0.999 

2.181 

0.984 

0.998 

2.184 

0.983 

0.999 

2.183 

0.985 

0.998 

2.186 

0.986 

0.998 

2.187 

0.985 

0.998 

2.186 

0.986 

0.998 

2.188 

0.988 

0.998 

2.189 

See  Table  5 for  crystal  aggregate  properties. 


Figure  11.  Concentration  ratio  vs.  particle  mass  for  various  hydrometeors  at  5 km  altitude. 

For  the  heavy  particle  asymptote,  CM  E « n_1  = 2.222.  Tables  2-6  list  particle 
masses  for  the  hydrometeors  studied.  Data  for  plane  dendrites  and  aggregates  are 
not  Included  since  they  fall  between  the  curves  shown  on  the  already  crowded  right- 
hand  side  of  the  graph. 


From  Fig.  11  and  Tables  6 - 10  we  see  that  the  largest  deviation  of 
concentration  ratio  from  the  heavy  particle  limit  Is  about  ten  percent  for 
particles  heavier  than  0.5  yg.  For  smaller  particles,  the  concentration 
ratio  falls  away  from  the  heavy  particle  limit  toward  unity,  with  the 
smallest  value  calculated  being  1.35  for  the  case  of  hollow  columns  of 
20  ym  length  at  2 km  altitude. 
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